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I.  INTRODUCTION 


It  has  been  known  for  quite  some  time  that  the  near-resonant  interaction 
between  electromagnetic  fields  and  atomic  or  molecular  systems  can  lead  to 
level  shifts  and  splittings  of  the  system  under  study. *  When  the  field  is 
sufficiently  weak,  the  shift  of  the  energy  levels  is  quadratic  in  the  field 
strength;  when  the  field  is  strong,  the  shift  is  linear  in  the  field 
strength.  The  correct  dependence  of  the  shift  on  field  strength  is  very 
important  when  one  performs  high-resolution  spectroscopy,  or  uses  the  field 
effect  to  accurately  measure  transition  dipole  matrix  elements. ^  However,  it 
is  equally  important  to  consider  the  effects  of  inhomogeneities  when 
performing  these  measurements,  since  optically  produced  inhomogeneous 
broadening  can  drastically  affect  the  observed  shift. 

Inhomogeneous  effects  can  be  easily  observed  and  studied  in  hyperfine 
optical  pumping  experiments  of  alkali  atoms,  where  the  atom  can  be  considered 

O 

as  essentially  "frozen  in  place"  because  of  the  presence  of  buffer  gas.  In 
these  experiments  the  level  shift  is  referred  to  as  the  light  shift  due  to 
virtual  transitions,  since  the  optical  field  produces  a  shift  observed  in  a 
microwave  or  rf  resonance.  These  shifts  were  first  discussed  by  Barrat  and 
Cohen-Tannoud j i , ^  and  then  treated  semiclassically  by  Happer  and  Mathur.-*  The 
light  shift  is  of  considerable  importance  in  a  number  of  devices  such  as  fre¬ 
quency  standards,  optically  pumped  magnetometers,  and  masers;  and  although 
inhomogeneous  effects  have  been  discussed no  comprehensive  experimental 
study  of  an  inhomogeneous  light  shift  has  been  undertaken. 

We  report  here  measurements  on  the  level  shifts  in  the  ground-state 
87 

Rb  (rubidium)  hyperfine  levels  induced  by  laser  radiation  approximately  in 
resonance  with  the  transition  from  the  ground  state  to  the  first  excited 
state,  5  Sj/2  “  5  Pj/2*  These  shifts  are  deduced  by  measuring  the  center 
frequency  of  the  ground-state  hyperfine  resonance  as  a  function  of  incident- 
laser  intensity.  When  inhomogeneous  broadening  is  present,  observed  as  an 
asymmetry  in  the  hyperfine  resonance  line,  the  center  frequency  is  found  to  be 


II.  EXPERIMENTAL  PROCEDURE 


The  experimental  apparatus  is  shown  schematically  in  Fig.  1.  A  glass 

absorption  cell  that  contains  an  excess  of  Rb  metal  and  ^  at  10  Torr  is  situ- 

ated  in  a  TEjjj  microwave  cavity  tuned  to  the  Rb  ground-state  hyperfine 
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transition  frequency.  The  Rb  is  in  its  natural  isotopic  abundance,  27.2 % 

®^Rb  and  72.8%  ®^Rb;  the  N2  is  present  in  order  to  quench  the  Rb  fluorescence, 
and  to  act  as  a  buffer  to  reduce  the  effect  of  collisions  with  the  cell  walls. 
A  static  magnetic  field  of  about  300  mG  is  applied  parallel  to  the  cavity  axis 
in  order  to  define  the  quantization  axis,  and  to  split  the  Zeeman  levels  so 
that  only  the  0-0  transition  is  induced  by  the  microwaves.  The  cavity  and 
cell  are  thermostatically  controlled  to  ±0.1°C  at  about  706C,  and  are  sur¬ 
rounded  by  a  single  layer  of  magnetic  shielding.  A  single  magnetic  shield  is 
sufficient,  since  the  transition  of  interest  is  independent  of  magnetic  field 
to  first  order. 

The  absorption  cell  is  illuminated  by  the  emission  from  a  single-mode 
diode  laser,  Mitsubishi  ML-4001 ,  that  is  tuned  to  one  of  the  87Rb  Dj  hyperfine 
resonance  lines  (794.7  nm);  this  produces  a  hyperfine  polarization  in  the 
ground  state.  A  microwave  field  is  either  swept  or  modulated  through  the  87Rb 
hyperfine  resonance  at  approximately  6834  MHz,  and  the  ensuing  change  in  the 
light  transmitted  through  the  absorption  cell  is  monitored  by  a  silicon  photo¬ 
diode. 

The  experimental  apparatus  is  configured  in  two  different  ways,  depending 
upon  whether  we  are  measuring  the  light  shift  or  the  line  shape.  These  two 
configurations  are  represented  in  Fig.  1  by  the  two  switch  positions  "A"  and 
"B."  These  are  not  physical  switches,  but  are  displayed  in  the  figure  in 
order  to  facilitate  the  understanding  of  our  set-up.  The  rationale  behind 
using  two  different  configurations  was  simply  experimental  convenience.  The 
voltage-controlled  crystal  oscillator  (VCX0)  frequencies  shown  in  Fig.  1  are 
not  crucial,  but  are  what  we  had  available  to  perform  the  different  measure¬ 


ments. 
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Experimental  Arrangement.  Switch  positions  A  and  B  correspond 
to  the  arrangements  used  to  measure  the  light  shift  and  line 
shape,  respectively. 


In  measuring  the  magnitude  of  the  light  shift,  one  needs  to  determine,  to 
a  high  degree  of  precision,  the  center  of  the  optically  detected  hyperfine 
line  as  a  function  of  the  intensity  and  spectral  detuning  of  the  pumping 
radiation.  To  perform  this  measurement,  we  frequency-lock  a  nominal  10-MHz 
VCXO  to  the  peak  of  the  hyperfine  transition  line.  This  is  schematically 
displayed  in  Fig.  1  when  the  switches  are  in  the  "A"  position.  The  output 
frequency  of  the  VCXO  is  multiplied  up  to  60  MHz  by  means  of  standard  cir- 
cuitry,  and  then  by  a  "step-recovery"  diode  to  6840  MR2  in  the  microwave 
cavity.  Since  this  is  not  the  proper  frequency  for  inducing  hyperfine  transi¬ 
tions,  a  synthesizer  is  used  to  provide  a  frequency  of  approximately  5  MHz, 
which  is  subtracted  from  the  6840  MHz  to  yield  the  correct  hyperfine  resonance 
frequency.  By  phase-modulating  the  microwave  radiation,  an  error  signal 
generated  at  the  photodiode  is  fed  back  to  the  VCXO.  This  scheme  is  similar 
to  the  typical  method  for  locking  a  VCXO  to  the  atomic  transition  in  a  Rb 
frequency  standard.*®  The  VCXO  frequency  is  averaged  for  100  s,  and  then 
read  with  a  frequency  counter.  From  the  circuitry,  the  relationship  between 
the  VCXO  frequency  and  the  microwave  frequency  is  known.  With  this  cali¬ 
bration  we  are  able  to  determine  the  hyperfine  resonance's  peak  frequency  to 
about  1  Hz. 

The  line  shape  of  the  optically  detected  hyperfine  transition  is  measured 
by  ramping  the  frequency  of  a  VCXO  through  the  hyperfine  resonance  and  moni¬ 
toring  the  light  transmitted  through  the  absorption  cell.  We  employ  a  linear 
VCXO  that  has  been  carefully  calibrated,  allowing  us  to  convert  VCXO  voltage 
directly  to  frequency.  The  microwave  field  is  chopped  with  a  diode  switch  at 
a  few  Hz,  and  a  lock-in  amplifier  is  used  to  enhance  the  signal-to-noise  ratio 
on  each  sweep.  Typically,  16  repetitive  scans,  taking  about  2  min  each,  are 
averaged  with  the  aid  of  an  HP  9825  computer.  We  estimate  an  accuracy  of 
about  5  Hz  in  the  line  shapes  with  this  method. 

We  do  not  have  a  direct  measure  of  the  strength  of  the  microwave  field  in 
the  cavity.  Multiplication  up  to  the  microwave  frequency  from  the  nominal 
VCXC  frequency  is  accomplished  using  a  "step-recovery"  diode.  Ideally,  the 
microwave  power  in  the  Nt*1  harmonic  is  expected  to  be^ 

PN  -  P0/N  (1) 
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where  Pq  is  the  input  power  to  the  step-recovery  diode.  However,  the  coupling 
coefficients  between  the  diode  and  the  cavity  antenna,  and  between  the  cavity 
antenna  and  the  cavity,  will  modify  Eq.  (1).  We  believe  that  a  reasonable 
estimate  of  the  microwave  power  in  the  cavity  is  on  the  order  of  10  pW. 

The  spectral  linewidth  of  the  ML  4001  diode  laser,  measured  with  a  Fabry- 
Perot  interferometer,  was  found  to  be  about  100  MHz.  The  laser's  wavelength 
could  be  tuned  to  either  of  the  Doppler-broadened  hyperfine  absorption  lines, 
AVp  ~500  MHz,  by  varying  either  the  diode  temperature  or  the  injection 
current.  The  diode  laser  is  heat-sunk  into  a  copper  block  whose  temperature 
is  stabilized  and  controlled  by  a  thermistor  in  one  leg  of  a  bridge  circuit 
that  controls  the  current  through  a  solid-state  thermoelectric  device.  In 
this  manner  the  diode  laser's  center  frequency  can  be  held  to  less  than  100 
MHz  of  the  center  of  the  hyperfine  absorption  line  for  about  30  min  without 
active  stabilization  of  the  laser  diode. 

Typically,  the  laser  is  tuned  by  first  adjusting  the  temperature  so  that 
the  lasing  wavelength  is  near  the  Rb  Dj  line  at  794.7  nm.  The  injection 
current  is  then  used  as  a  fine  control  to  tune  the  laser  over  the  hyperfine 
absorption  spectrum,  and  as  long  as  lasing  mode  hops  do  not  occur,  the  injec¬ 
tion  current  can  be  calibrated  to  the  lasing  frequency.  This  is  found  to  be 
approximately  16  GHz/mA.  Since  the  laser's  single-mode  output  power  is  a 
function  of  the  injection  current,  the  laser  power  will  vary  as  the  laser  is 
tuned.  However,  the  fractional  change  in  power  is  found  to  be  only  about  1% 
as  the  laser  is  tuned  over  several  GHz,  and  thus  the  variations  in  laser  power 
are  neglected  for  these  measurements.  The  typical  total  laser  power  in  the 
single-mode  line  is  found  to  be  about  3  mW  in  a  Gaussian  beam  diameter  of 
MJ.45  cm. 
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III.  RESULTS 


The  transmission  of  the  laser  diode  emission  as  a  function  of  the  laser's 
frequency  is  shown  in  Fig.  2,  along  with  a  labeling  of  the  observed  transi¬ 
tions.  Doppler  broadening  prevented  resolution  of  the  excited-state  hyperfine 
structure,  except  for  the  ®^Rb  transitions:  *  O  “  5Pj /2 ^ F  2). 

By  pumping  the  high-frequency  component  of  these  transitions,  the  effects  of 
®^Rb  absorption  were  minimized,  and  except  where  noted  all  detunings  are 
referenced  with  respect  to  this  component. 

A.  LIGHT-SHIFT  MEASUREMENTS 

The  light  shift  can  be  interpreted  as  being  due  to  either  non-energy- 
conserving  virtual  transitions,11  or  as  the  interaction  between  the  induced 
polarization  of  the  atom  and  the  electric  field  of  the  light. ^  Both  of  these 
interpretations  are  consistent  with  second-order  perturbation  theory,  which 
predicts  a  linear  relationship  between  the  shift  and  the  light  intensity.  In 
either  case,  one  can  show  that  the  light  shift  can  be  expressed  qualitatively 
as 


Av 


hf  s 


(2) 


where  constants  have  been  omitted  for  clarity.  The  plus  and  minus  signs  refer 
to  excitation  from  the  lower  or  upper  ground-state  hyperfine  level,  respec¬ 
tively.  The  incident-light  intensity  is  denoted  by  Iq,  the  atomic  transition 
frequency  is  denoted  by  Vq,  the  laser  emission  frequency  is  v^»  and  t  is  the 
spontaneous  lifetime  of  the  excited  state. 

Our  results,  shown  in  Fig.  3,  of  the  light  shift  versus  incident-laser 
intensity  are  markedly  different  from  the  strict  linear  dependence  predicted 
by  the  above  expression.  From  top  to  bottom,  the  four  curves  represent  laser 
detunings  of  approximately  -400  MHz,  -200  MHz,  0,  and  +200  MHz.  All  of  these 
curves  show  the  same  general  behavior:  they  are  linear  at  low  laser  inten¬ 
sity,  reach  an  extremum,  and  finally  saturate  to  a  light-shift  value  that 
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Rb  Hyperfine  Transitions  Observed  in  Transmission  of  the  Laser 
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seems  to  depend  only  on  the  laser  detuning  from  the  atomic  resonance.  The 
slight  light  shift  observed  when  the  laser  is  tuned  on  resonance  can  be  easily 
explained  by  the  partial  overlap  of  the  5Sj^^  “  D  ~  5P ^ / 2 ^ F  =  1»  2) 
transitions. 

The  linearity  of  these  curves  at  low  light  intensity,  and  the  sign  of 
their  symmetry  about  zero-frequency  detuning,  are  consistent  with  the  above 
expression.  However,  Eq.  (2)  does  not  predict  either  the  extrema  or  the 
saturation  regions  in  the  light-shift  curves.  The  extrema  and  saturation 
regions  are  also  not  consistent  with  a  dressed  atom  model  of  the  light  shift 
which  predicts  a  square-root  dependence  of  the  shift  on  high  light  inten- 

1  -J  O 

sity,  as  demonstrated  by  Liao  and  Bjorkholm.  In  order  to  understand  this 
discrepancy  more  fully,  we  undertook  a  second  set  of  experiments  to  analyze 
the  line  shape  of  the  microwave  resonances  as  a  function  of  laser  intensity 
and  tuning. 

B.  LINE  SHAPE  MEASUREMENTS 

In  Fig.  4  we  present  representative  samples  of  our  microwave  line  shape 

measurements.  The  upper  line  shapes  correspond  to  low  incident-laser  inten- 
2 

sity,  ~0.5  mW/cm  ,  for  several  different  detunings  of  the  laser  frequency. 

From  left  to  right  these  are  -A20,  0,  and  +420  MHz.  The  lower  curves  were 
observed  with  full  laser  intensity,  ~10  mW/cm  .  From  the  figure  it  is  ap¬ 
parent  that  for  higher  laser  intensities  the  microwave  line  shape  is  asym¬ 
metric  when  the  laser  is  tuned  off  resonance,  and  that  this  asymmetry  has  the 
same  sign  as  the  light  shift. 

To  summarize,  then,  we  find  that  with  high  laser  power  the  light  shift 
has  a  nonlinear  dependence  on  light  intensity,  and  that  this  nonlinearity 
appears  to  be  connected  with  an  optically  induced  asymmetry  In  the  microwave 
line  shape.  From  these  observations  one  might  be  tempted  to  assume  that  these 
effects  are  due  to  the  onset  of  dynamic  Stark  splitting.  (In  our  experiment 
nr  ~  0.1,  where  SI  is  the  optical  Rabi  frequency.)  However,  the  nonlinearity 
does  not  have  the  expected  square-root  dependence,  and  the  sign  of  the 
microwave  asymmetry  is  wrong:  for  an  unresolved  splitting  one  would  predict 
two  asymmetric  resonances  shifting  away  from  one  another  with  increasing  laser 
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5  mW/cm 


Experimental  Result*  of  the  Microwave  Line  Shape 
Measurements.  The  top^row  Is  for  low  laser 
Intensity,  ~  0.5  mW/cm  ,  and  the  bottom  row  for  high  laser 
Intensity,  ~  10  mW/cm^.  From  left 

to  rtp.ht  the  laser  detuning*  are  -470,  0,  and  +470  M07.. 
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intensity,  so  that  the  sign  of  the  asymmetry  would  be  opposite  that  of  the 
observed  light  shift.  In  the  following  section  we  show  that  these  effects  can 
be  adequately  described  by  the  inhomogeneous  light  shift,  which  does  not 
necessarily  include  the  phenomenon  of  dynamic  Stark  splitting. 


IV.  MODEL  OF  THE  INHOMOGENEOUS  LIGHT  SHIFT 


In  general,  different  atoms  in  a  gas  experience  different  local  perturba¬ 
tions  of  their  atomic  states.  For  the  hyperfine  states  of  ®^Rb,  as  shown  in 
Fig.  %  the  differences  are  usually  the  result  of  gradients  over  the  cell 
dimensions:  temperature  gradients,  static  magnetic-field  gradients,  microwave 

energy-density  gradients,  light-intensity  gradients,  etc.  If  these  gradients 
are  large  enough,  and  if  the  atom  cannot  average  over  them,  then  the  hyperfine 
transition  will  be  inhomogeneously  broadened. 

In  an  attempt  to  model  these  inhomogeneous  effects,  consider  a  cylindri¬ 
cal  cell  divided  into  a  number  of  finely  spaced  layers,  where  protruding 
through  each  layer  are  a  number  of  finely  spaced  volume  elements.  We  assume 
that  the  atoms  within  a  volume  element  all  experience  the  same  perturbations, 
and  can  be  described  by  a  local  density  matrix  determined  by  rate  equations 
similar  to  those  used  by  Missout  and  Vanier:^ 

p(r ,6  ,z)  =  p(r,0,z)relax  +  p(r,0  ,z)rf  +  p(r,0 ,z)opt  (3) 

The  three  terms  on  the  right-hand  side  of  Eq.  (3)  correspond  to  the  effects  of 
relaxation,  the  microwave  field,  and  optical  excitation,  respectively;  each  of 
these  will  be  discussed  more  fully  below.  In  addition,  we  assume  that  the 
microwave  transition  is  Dicke  narrowed,^  so  that  Doppler  effects  in  the 
microwave  spectra  can  be  ignored. 

The  local  density  matrix  determines  the  fraction  of  light  transmitted  by 
a  particular  volume  element.  Rows  of  these  volume  elements  determine  the 
fraction  of  incident  light  transmitted  by  the  cell.  We  perform  the  calcula¬ 
tion  by  first  considering  a  particular  row.  We  assume  we  know  the  light 
intensity  incident  on  this  row,  and  average  it  over  the  row's  cross-sectional 
area.  We  also  assume  that  we  know  the  microwave  magnetic  field  strength  at 
any  point  along  the  row,  and  average  it  over  a  given  volume  element.  We  can 
then  step  through  the  row  of  volume  elements,  calculating  the  average 


F  =  2 


light-intensity  attenuation  for  each  step.  In  this  way  we  can  determine  the 
light  intensity  transmitted  through  the  cell  as  a  function  of  microwave 
frequency,  thus  obtaining  the  microwave  line  shape. 


Once  the  line  shape  has  been  calculated,  it  is  a  relatively  simple  matter 
to  determine  the  frequency  corresponding  to  the  peak  of  the  line  shape;  this 
is  the  frequency  that  is  measured  in  our  experiments.  We  should  note  that 
this  is  equivalent  to  determining  the  resonant  frequency  for  the  homogeneous 
packet  of  atoms  that  has  the  greatest  contribution  to  the  line  shape.  When 
local  perturbations  differ  significantly  from  one  another,  we  refer  to  the 
dependence  of  this  resonant  frequency  on  the  incident-light  intensity  as  the 
inhomogeneous  light  shift,  in  order  to  distinguish  it  from  the  more 
fundamental  light  shift  discussed  by  Happer  and  Mathur.^ 

A.  THE  LOCAL  DENSITY  MATRIX 

The  first  term  in  Eq.  (3)  refers  to  uniform  relaxation  of  the  density 
matrix.  Since  our  present  interest  is  in  the  effects  of  light-intensity 
gradients,  we  assume  that  the  rate  constants  for  this  term  are  independent  of 
spatial  position  within  the  cell.  Following  Missout  and  Vanier,*-*  we  express 
this  term  as 


iiw<r.e.z>]relax  ■  -Vi  -  <,vv(r'e,z)l  <*> 

[;vv,(r,e,z)]relai(  -  Kj  Pvv,(r,e,.)  <v*v’)  (5) 

where  we  have  assumed  that  only  two  rates  are  necessary  to  describe  the 
relaxation  of  the  ground-state  multiplets  and  yj,  which  refer  to  the  long¬ 
itudinal  and  transverse  relaxation  rates,  respectively. 

The  second  term  of  Eq.  (3)  describes  the  coherences  produced  in  the 
density  matrix  as  a  result  of  the  application  of  microwaves  at  the  hyperfine 
transition  frequency.  Since  the  cell  is  situated  in  a  TEj^  microwave  cavity, 
this  term  will  have  a  very  sensitive  dependence  on  spatial  position.  We  will 
simplify  our  calculation  by  assuming  that  the  only  coherence  generated  in  the 


density  matrix  is  between  the  (F,mp)  (2,0)  -  (1,0)  states.  Since  the  Zeeman 
splitting  is  ~100  times  larger  than  the  width  of  the  observed  line  shapes,  oar 
assumption  corresponds  quite  closely  to  our  experimental  situation.  Thus,  we 
only  need  to  consider  the  terms 


[p77(r,e  ,z)  Jrf  -  2ojj  (r  ,6  ,z)  1m  [p  37(r  ,6  ,z)eiui  C] 


[p33(r,6  ,z)]rf  -  -2u)j  (r ,6  ,  z)  Im  [p  3?(r  ,0  .zje1^] 


[p3?(r,e,z)]rf  =  io)1(r,e,z)[p33(r,6,z)  -  p  7?(r,0  ,z)]  e" 


-  iajQp37(r ,0,z) 


where  is  the  unperturbed  hyperfine  resonance  frequency,  and  u j(r,6  ,z)  is 
the  spatially  dependent  microwave  Rabi  frequency. 

The  last  term  in  Eq.  (3)  corresponds  to  the  optical  pumping  process. 

This  term  has  a  radial  dependence  that  results  from  the  intensity  profile  of 
the  laser  beam,  and  an  axial  and  angular  dependence  that  results  from  optical 
absorption  by  successive  layers  of  Rb  atoms  in  the  cell.  Since  the  laser 
linewidth  is  much  narrower  than  the  hyperfine  splitting,  but  much  broader  than 
the  Zeeman  splitting,  the  laser  pumps  atoms  out  of  only  one  hyperfine  state, 
but  out  of  all  Zeeman  sublevels  of  that  state  equally.  As  shown  in  Fig.  3,  we 
consider  the  optical  transitions  5S.  ,_(F  -  1)  -  5P.  ,«(F  ■  1,  2),  so  that  from 
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Missout  and  Vanier1 J  we  have 


[p  (r,6,z)]  _  ■  -r(r,0,z)p  (r,0,z)6 

L  vv  Jopt  vv  pv 


♦  !  p  <r,e ,,)  (,-6.7.8) 
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(H) 


•c37(r,£,z) -opt 


r  r(r ,5 ,z) 
2 


i6w(r ,6  ,z) ^  p37(r,6  ,z) 


where  r(r,0,z)  and  6w(r,6,z)  are  the  spatially  dependent  pumping  rate  and 
light  shift,  respectively. 

Two  important  assumptions  are  made  in  writing  Eqs.  (9)  and  (10).  First, 
because  of  the  presence  of  a  significant  amount  of  nitrogen  in  our  cell,  we 
assume  that  relaxation  from  the  P  state  is  predominantly  nonradiative , ^  and 
that  it  occurs  with  equal  probability  to  all  of  the  ground-state  sublevels. 
Second,  since  the  optical  Rabi  frequencies  that  we  achieve  in  our  experiment 
are  much  smaller  than  the  relaxation  rate  from  the  excited  P  state,  we  assume 
there  is  negligible  mixing  of  the  ground  and  excited  states.  Thus,  this  model 
of  the  inhomogeneous  line  shape  is  essentially  a  low-light-intensity  model. 

Since  we  are  only  concerned  with  relatively  low  light  intensities,  we 
expect  the  light  shift  to  be  proportional  to  the  light  intensity  at  any  point 
in  the  cell,  as  predicted  by  second-order  perturbation  theory.  This  is  quite 
convenient  for  our  computation,  since  we  need  only  calculate,  using  the  theory 
of  Happer  and  Mathur,^'^  the  light  shift  once.  However,  this  will  only  be 
correct  if  the  spectral  line  shape  of  the  laser  does  not  change  as  the  beam 
propagates  through  the  cell.'**-’  In  the  present  experiment,  the  optical 
absorption  line  is  Doppler  broadened  to  ~500  MHz,  so  that  we  wouldn't  expect 
the  ~100  MHz  laser  line  to  experience  any  appreciable  distortion  as  it  passes 
througn  the  vapor. 

From  the  symmetry  of  our  model  it  is  apparent  that 


P11  "  p22  “  p44  "  p55  ;  P66  “  p88 


(ID 


Thus,  we  need  only  consider  solutions  to  the  following  equations: 

Pn(r,9,z)  [p  (r,8,z)+71[J-p11(r,0,i)]  (12) 

v 
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P33(r,e,z)  =  ’g’  z  pyi/r»9  *z)  +  Y  ll-  8  "  p33(r,e  ,z)j 

P 

-  ZidjCr.e.z)  Im  [p3?(r,e  ,z)elwt] 


p66(r,e,z)  =  -r(r,e,z)  p66(r,e,z)  +  r(r,g,z)  E  Pyy(r,e,z)  + 


Vf"  p66(r’6'z)^ 


P77(r,0,z)  ■=  -r(r,6,z)  p77(r,0,z)  +  ISljfLill  z  Pyy(r,8,z)  + 


T  i  [■g"  —  P77(r,0,z)]  +  2uij(r,0,z)  Im  [p  37(r,0  .zje1^1']  (15) 


p37(r,0,z)  -  -[(r(r,7,z)  +  Y2)  +  +  «u(r,0,z))]  P37(r,0  ,z)  + 


ioijCr.O.z)  [p33(r,0,z)  -  P7?(r,0  .z^e-*"*  (16) 

where  y  *  6,  7,  8. 

If  we  assume  a  solution  for  p37(r,0,z)  of  the  form  p  37(r,0  ,z)e  , 
where  p37(r,0,z)  may  be  a  slowly  varying  function  of  time,  then  in  steady- 
state  we  have 

Im  [p37(r,8,z)eiu,t]  -u^r.O.z)  c(r,0,z)  [p33(r,0,z)  -p7?(r,0,z)]  (17) 


where  we  call  c(r,0,z)  the  local  line  shape  factor;  i.e., 

[rir^  +  Ti] 


Using  Eq.  (17)  in  Eqs .  (12)  through  (15)  we  can,  after  some  algebra, 
obtain  analytical  expressions  for  the  local  steady-state  density  matrix 
elements  : 


P77(r,6,z)  = 


Yj[Yj  +  4tu1  (r,e,z)  c(  r ,  8,z  )  ] 


A(r , 6,z) 


Y1r(r,e,z) 

P33(r,e,z)  =  p77(r,e,z)  +  ■■  X(r,e7z-) 
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a  _x  1  rr(r’e,z)  p77  (r’e*z)  +  *1, 

(r’  ,z)  ■  9  [  3r(r,e,z)  +  4y  J 


1 


by, 


where 
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Pii(r,e,z)  =  i  +  [?■%£&>]  p66(r,e,z)  +  [•"(r-8^;z)]  P77(r,e,z)  (22) 
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A(r ,  6,z )  -  4iDj2(r,6,z)  c(r,6,z)  [r(r,6,z)  +  8Yj  +  e^zj’+^Y  i  + 


Y1[5r(r,e,z)  +  8Yj] 


(23) 


B.  THE  MICROWAVE  RABI  FREQUENCY 

In  order  to  determine  explicitly  the  local  density  matrix  elements,  it  is 
necessary  to  have  an  expression  for  the  local  Rabi  frequency.  Since  we  are 
only  considering  the  0-0  hyperfine  transition,  the  Rabi  frequency  is  given  by 


u^r.e.z) 


(r,6,z) 


(24) 


where  is  the  Bohr  magneton  and  Bz  is  the  axial  component  of  the  microwave 
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magnetic  field.  For  a  TEni  cavity,  this  field  component  has  the  form 


(25) 


Bz(r,6,z) 
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cos6 
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sin,— 


where  R  and  L  are  the  cavity  radius  and  length,  respectively.  The  field 
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amplitude  Bq  can  be  related  to  the  energy  density  in  the  cavity,  and  thus  to 
the  cavity  Q  and  power  loss.  After  performing  the  necessary  algebra,  we  have 
for  the  local  Rabi  frequency 


1/2 


4 .<r,9,z)  .12°  [O  ->2/4L2)_4S.<P>_]-,»  j  ( i^li)  cos6  si<^  (26) 

1  tt  R^  L  (0.236)v0  1 


where  <P>  is  the  rms  microwave  power  supplied  to  the  cavity. 


Averaging  Eq.  (26)  over  the  volume  element  r6r666z,  we  have,  for  the 
average  local  Rabi  frequency  used  in  the  calculation,  the  following: 


<(i)j(r ,0,z)>  -  tc  [ 


4L  sin(68/2)  6in(^-|^) 


■]  cos6  sin[^-)  x 
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and 


An  “  (_1)Tl  (~^-)2n+l  [2Zn+1(2n+3)  n!  (n+1)!]”1  (29) 

C.  THE  PUMPING  RATE  AND  LIGHT  SHIFT 

Both  the  puaping  rate  and  light  shift  are  functions  of  the  light  inten¬ 
sity  and  spectral  overlap  of  the  laser  line  and  Doppler-broadened  absorption 
line.  In  particular,  we  have  for  the  local  pumping  rate 
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r(r  ,8  ,z)  =  /  I(r  ,6  ,z,v  )  oD  (v  )  d*; 

where  v  is  the  center  frequency  of  the  laser.  For  convenience,  we  will 
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assume  that  the  laser  line  is  Gaussian,  7  so  that 

.  —  v  U  ftn2i  1/2  r  ..  ~  r  ’^L'l  2i 

I(r  ,0  ,z,v)  =  I(r,8,z)  —  [— ]  exp[-  4in2  (-^-j  j 
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where  Av  is  the  laser  linewidth  and  l(r,0,z)  is  the  total  local  intensity. 
Using  Eq.  (31)  in  (30)  we  have,  for  the  local  puraping  rate, 


r(r,e,z)  =  I(r,e,z)411  —  [ — 1T-/n2 - j]  exp[-  Ui.nl  ( 

0  Avl  +  Avd 


2  2 
Avl  +Avd 


]  (32) 


where  Vq  is  the  center  frequency  of  the  optical  absorption  line;  AvQ  is  the 
absorption  linewidth;  5v  is  the  detuning  of  the  laser  from  the  absorption 
line;  and  Tq,  f,  and  c  are,  respectively,  the  classical  electron  radius,  the 
oscillator  strength  of  the  transition,  and  the  speed  of  light. 

From  our  measurements  we  can  approximate  the  laser  with  a  Gaussian  radial 
intensity  distribution  at  the  entrance  of  the  cell: 


I(r,6,0)  =  —  PQ  exp[  -  4tn2  (r/a)2] 


where  a  is  the  Gaussian  beam  diameter  and  Pq  is  the  total  power  in  the  beam. 
We  average  Eq.  (33)  over  an  arbitrary  row’s  crosB-aectional  area,  which  gives 


<l(r,6,0)> 


PQ  exp[-  4Jtn2(r/a)2] [  1  -  ln2  (6r/a)2] 


As  the  beam  propagates  through  a  row  in  the  cell,  each  volume  element 
attenuates  the  beam  by  an  amount  <AI(r,6 ,z)>^.  Thus,  the  average  local 
intensity  at  the  nc^  volume  element  of  a  particular  row  is 


(35) 


<I(r,e,z)>  =  <I(r,6,0»  -  Z  <Al(r,6 ,z)>. 

i*l 

If  we  assume  exponential  attenuation  of  the  light  as  it  passes  through  the 
vapor,  then 


<AI(r,e,z)>i  -  -  n1  [Rb]  <r(r,e,z»i  2^62  (361 

where  is  the  average  fraction  of  atoms  in  the  i**1  volume  element  which  can 


absorb  light,  so  that 

ni  *  <p66^r,e,Z^>i  +  <P77(r»e.z)>i  +  <Pg8(r,e,z»i  (37) 

Equation  (35)  can  thus  be  used  in  Eq.  (32)  to  compute  the  average  local 
pumping  rate  in  the  nc^  volume  element,  and  also  the  average  local  light  shift 
for  the  volume  element: 


<6w(r,6 ,z»n  -  6*0  [<I(r,e,z)>n 


where  6u>n  is  the  light  shift  corresponding  to  a  light  intensity  of 


PQ  HL  n2 


V.  DISCUSSION 


In  order  to  obtain  a  clear  picture  of  the  inhomogeneous  light  shift,  we 
have  performed  two  sets  of  calculations.  In  the  first,  we  consider  the  in¬ 
homogeneous  effects  in  a  TEjjj  microwave  cavity;  this  should  correspond 
closely  to  our  experimental  data.  In  the  second  set,  we  imagine  a  fictitious 
cavity  where  the  microwave  magnetic  field  strength  is  everywhere  constant,  and 
directed  parallel  to  the  cavity  axis.  Microwave  line  shapes  and  the  light 
shift  versus  input  laser  intensity  for  both  microwave  field  configurations  are 
shown  in  Figs.  6  through  9.  The  input  parameters  for  these  calculations  are 
listed  in  Table  1. 

In  Fig.  6  we  compare  the  microwave  line  shapes  calculated  for  a  TE^^ 
microwave  cavity  at  low  and  high  laser  intensities.  When  the  laser  intensity 
is  low,  the  line  shape  appears  quite  symmetric.  This  is  because  the 
conditions  for  light-induced  inhomogeneous  broadening,  as  discussed  in  the 
Appendix,  are  not  satisfied  at  very  low  laser  intensity.  When  the  laser 
intensity  is  high,  the  line  shape  is  highly  asymmetric,  with  the  asymmetry  in 
the  direction  of  the  light  shift. 

The  sign  of  the  asymmetry  is  related  to  the  transmission  signal  amplitude 

for  the  local  volume  elements.  In  the  limit  of  high  light  intensity 

2  20 

(T  »  ojj),  this  amplitude  is  proportional  to  (u^/D  .  Thus,  the  homogene¬ 
ous  packets  of  atoms  that  experience  the  largest  light  shift  contribute  the 
least  to  the  total  signal  amplitude.  In  the  TEjn  microwave  cavity,  this 
effect  is  further  enhanced  by  the  overlap  of  the  radial  components  of  the 
microwave  and  optical  fields.  The  optical  field  has  a  radial  distribution 
described  by  a  Gaussian  function,  the  microwave  field  a  radial  distribution 
described  by  a  first-order  Bessel  function.  Thus,  for  radial  positions  where 
the  optical  field  is  most  intense,  the  microwave  field  is  not,  and  vice  versa. 

In  Fig.  7  we  compare  the  low-  and  high-intensity  line  shapes  when  the 
microwave  field  is  constant.  Again,  for  low  laser  intensity  the  line  shape  is 
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Calculated  Resonant  Microwave  Line  Shapes  in  a  TE^jj  Microwave 
Cavity  for  Low  and  High  Laser  Intensity* 
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Fig.  7.  Calculated  Resonant  Microwave  Line  Shapes  In  a  Fictitious 

Microwave  Cavity  with  Constant  rf  Magnetic  Field  for  Low  and 
High  Laser  Intensity. 
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INHOMOGENEOUS 
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8.  Calculated  Inhomogeneous  Light  Shift  In  a  TEj, j  Microwave 

Cavity  for  Several  Detunings  of  the  Laser.  From  top  to  bottom 
these  are  -411,  -201,  and  •♦■219  MHz. 
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AVERAGE  INPUT  LIGHT  INTENSITY  lmW/cmzl 


Calculated  Inhomogeneous  Light  Shift  in  a  Fictitious  Microwave 
Cavity  with  Constant  rf  Magnetic  Field  for  Several  Detunings  of 
the  Laser.  From  top  to  bottom  these  are  -411,  -201,  and 
+219  MHz. 


Table  1.  Calculation  Parameters 


laser  linewidth 
laser  beam  diameter 
absorption  linewidth 
longitudinal  relaxation  rate 
transverse  relaxation  rate 
rubidium  density 
cell  length 
cell  radius 

cavity  aperture  diameter 
cavity  quality  factor 
rms  microwave  power 


«  100 

MHz 

0.45 

cm 

-  500 

MHz 

100 

Hz 

o 

o 

Hz 

b]  - 

1  x  10 

3.8  cm 

1.35 

cm 

1.6  cm 

Note:  In  this  cavity,  10  yW  of  rf  input  power  implies  an  average  axial 
magnetic  field  of  approximately  9  yG. 


symmetric,  and  at  high  laser  intensity  the  line  shape  is  asymmetric.  The  sign 
of  the  asymmetry  is  also  the  same  as  for  a  TE^j  cavity,  which  we  would  expect 
from  the  arguments  presented  above. 

In  the  Appendix,  the  necessary  conditions  for  observing  inhomogeneous 
broadening  are  discussed  in  terms  of  the  dephasing  rate  ( 1  /T2 )  an<^  the  longi¬ 
tudinal  relaxation  rate  (1/Tj).  However,  each  of  these  rates  is  a  function  of 
both  the  pumping  rate  and  the  microwave  Rabi  frequency.  Therefore,  each  of 
these  rates  varies  spatially,  and  as  a  consequence  so  do  the  degrees  to  which 
the  inhomogeneous  broadening  conditions  are  satisfied.  Thus,  the  inhomo¬ 
geneous  line  shape  is  not  composed  of  homogeneous  line  shapes  with  equal 
linewidths,  but  of  homogeneous  line  shapes  with  different  linewidths,  the 
linewidths  being  determined  by  the  local  T2's  and  the  degree  to  which  motional 
narrowing  occurs  on  a  local  level. 

Figure  8  shows  the  calculated  inhomogeneous  light  shift  as  a  function  of 
incident  light  intensity  in  a  TE ^ 1 1  cavity.  The  several  curves  correspond  to 
different  laser  detunings,  and  it  is  apparent  that  the  agreement  between 
theory  and  experiment,  shown  in  Fig.  3,  is  rather  good.  In  Fig.  9  we  show  the 
Inhomogeneous  light  shift  for  a  constant  microwave  field.  The  fact  that  the 
nonlinearity  is  observed  with  both  calculations  implies  that  this  feature  is 
primarily  determined  by  the  light-intensity  gradients  in  the  cell,  and  not  by 
the  microwave  field  distribution.  The  puzzling  thing  about  Figs.  8  and  9  is 
that,  though  extrema  in  the  inhomogeneous  light-shift  curves  were  observed 
experimentally,  such  extrema  are  only  predicted  in  the  constant-microwave- 
field  calculations. 

The  extrema  in  the  experimental  light-shift  curves  could  be  due  to  a 
number  of  effects  not  considered  in  our  calculation.  For  example,  the  0-0 
hyperfine  transition  frequency  has  a  quadratic  dependence  on  the  static 
magnetic  field  strength.  Thus,  magnetic  field  inhomogeneities  could  easily 
have  a  significant  effect  on  the  Inhomogeneous  light  shift.  However,  a  more 
intriguing  possibility  is  that  the  atoms  actually  do  experience  a  relatively 
constant  microwave  field.  This  might  occur  if  motional  narrowing  on  a  local 
level  allowed  some  averaging  of  the  microwave  and  optical  fields. 
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As  pointed  out  in  the  Appendix,  the  criterion  for  atoms  being  frozen  in 
place  might  not  be  strictly  fulfilled  everywhere  in  the  cell.  In  such  a  case 
the  microwave  and  optical  fields  might  be  averaged  in  some  regions  of  the  cell 
but  not  in  others.  The  results  of  this  could  be  to  effectively  change  the 
spatial  distribution  of  the  microwave  field  as  experienced  by  an  atom,  so  that 
the  "effective"  field  might  have  a  spatial  distribution  quite  different  from 
that  of  the  TE^jj  cavity.  If  this  local  motional  narrowing  is  occurring,  the 
extrema  in  the  light-shift  curves  should  be  a  function  of  the  buffer  gas 
density:  as  the  buffer  gas  density  is  increased,  the  diffusion  length  of  the 
atom  decreases,  and  less  averaging  of  the  fields  is  possible.  The  effective 
microwave  field  distribution  will  then  approach  the  TE^j  distribution  as 
the  buffer  density  is  increased,  and  the  measured  light-shift  curves  should 
approach  those  displayed  in  Fig.  8.  Unfortunately,  at  the  present  time  we  are 
unable  to  make  any  definite  statements  as  to  the  source  of  the  extremum. 


VI.  CONCLUSIONS 


We  have  found  that  inhomogeneous  broadening  of  a  spectral  line  as  a 
result  of  light-intensity  gradients  can  lead  to  a  nonlinear  dependence  of  the 
line  shape's  center  frequency  on  incident-light  intensity.  Thus,  one  needs  to 
be  very  careful  in  analyzing  optically  induced  energy-level  shifts,  either  for 
high-precision  spectroscopy  or  for  a  determination  of  transition  dipole  matrix 
elements,  since  this  nonlinearity  could  significantly  influence  the  results. 

In  addition,  one  should  be  very  careful  in  interpreting  asymmetries  in  reso¬ 
nance  lines  as  being  due  to  the  onset  of  an  Autler-Townes-type  splitting.* 
Except  for  the  sign,  asymmetries  due  to  Inhomogeneous  broadening  could  be 
easily  confused  with  this  effect.  It  is  our  contention  that  this  was  exactly 
the  case  with  the  asymmetries  observed  by  Arditi  and  Picque  for  the  0-0 
hyperfine  transition  of  cesium  optically  pumped  by  a  laser  diode.  Since  the 
intensities  that  they  could  have  reasonably  achieved  with  their  laser  diode 
were  probably  similar  to  ours,  we  believe  that  they  too  were  observing  the 
effects  of  inhomogeneous  broadening  as  discussed  in  this  report. 
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APPENDIX:  THE  POTENTIAL  FOR  MOTIONAL  AVERAGING  OF 
SPATIAL  INHOMOGENEITIES 


As  mentioned  in  the  text,  there  are  essentially  two  conditions  that  roust 
be  met  in  order  to  have  an  inhoroogeneously  broadened  microwave  line.  First, 
the  local  perturbation  variations  must  be  of  sufficient  magnitude,  so  that  the 
spread  in  perturbed  resonant  frequencies  is  greater  than  the  homogeneous 
linewidth : 


(6V£oc/h)  >  (1/T2}  (A1) 

where  represents  the  variation  in  local  perturbations,  and  I/T2  is  the 

homogeneous  linewidth.  In  addition,  motional  narrowing  should  be  insignifi¬ 
cant:  the  individual  atoms  should  not  be  able  to  average  the  local  perturba¬ 
tions.  In  the  present  case,  this  last  requirement  would  be  satisfied  if  the 
atoms  were  not  able  to  diffuse  spatially  over  the  effective  volume  of  local 
perturbation  gradients  during  the  lifetime  of  the  hyperfine  states. 

If  we  consider  diffusion  in  one  dimension,  then  during  the  time  Tj  an 
atom  will  diffuse  a  length 
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where  D  is  the  diffusion  coefficient  for  Rb  atoms  in  a  buffer  gas.  In  three 
dimensions,  we  can  think  of  the  atom  as  diffusing  through  a  volume 
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Thus,  for  motional  narrowing  not  to  occur  we  would  want 
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where  a  is  the  light  beam  radius  and  L  is  the  cell  length.  In  the  present 
case,  if  we  assume  a  lifetime  on  the  order  of  1  msec, 


Vdif/Veff  ~  5  *  10 


which  suggests  that  to  a  first  approximation  the  Rb  atoms  can  be  considered  as 
frozen  in  place  as  a  result  of  the  presence  of  10  Torr  of  N£.  Caution  must  be 
exercised  here  because  if  Tj  and  T2  vary  through  the  cell  volume,  then  the 
criterion  for  treating  the  atoms  as  frozen  in  place  might  not  be  fulfilled 
everywhere  in  the  cell.  In  that  case  some  local  averaging  of  the  perturba¬ 
tions  might  take  place  (i.e.,  local  motional  narrowing). 
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Computer  Science  Laboratory:  Program  verification,  program  translation, 
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computers,  fault-tolerant  computer  systems,  artificial  intelligence  and 
microelectronics  applications. 

Electronics  Research  Laboratory:  Microelectronics,  GaAs  low  noise  and 
power  devices,  semiconductor  lasers,  electromagnetic  and  optical  propagation 
phenomena,  quantum  electronics,  laser  communications,  Ildar,  and  electro- 
optics;  communication  sciences,  applied  electronics,  semiconductor  crystal  and 
device  physics,  radiometric  imaging;  millimeter  wave,  microwave  technology, 
and  RF  systems  research. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metal 
matrix  composites,  polymers,  and  new  forms  of  carbon;  nondestructive  evalua¬ 
tion,  component  failure  analysis  and  reliability;  fracture  mechanics  and 
stress  corrosion;  analysis  and  evaluation  of  materials  at  cryogenic  and 
elevated  temperatures  as  well  as  in  space  and  enemy-induced  environments. 

Space  Sciences  Laboratory:  Magnetospheric,  auroral  and  cosmic  ray  phys¬ 
ics,  wave-particle  Interactions,  magnetospheric  plasma  waves;  atmospheric  and 
ionospheric  physics,  density  and  composition  of  the  upper  atmosphere,  remote 
sensing  using  atmospheric  radiation;  solar  physics,  Infrared  astronomy, 
Infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
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